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We compute the Gilbert damping in (Ga,Mn)As based on the scattering theory of magnetization 
relaxation. The disorder scattering is included non-perturbatively. In the clean limit, spin-pumping 
from the localized d-electrons to the itinerant holes dominates the relaxation processes. In the 
diffusive regime, the breathing Fermi-surface effect is balanced by the effects of interband scattering, 
which cause the Gilbert damping constant to saturate at around 0.005. In small samples, the system 
shape induces a large anisotropy in the Gilbert damping. 



I. INTRODUCTION 

The magnetization dynamics of a ferromagnet can be 
described phenomenologically by the Landau-Lifshitz- 
Gilbert (LLG) equation:^ ^ 
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Here, 7 is the gyromagnetic ratio, Hcff is the effective 
magnetic field (which is the functional derivative of the 
free energy Hoff = —5F [M] /(5M), M is the magnetiza- 
tion and Mg is its magnitude. The Gilbert damping con- 
stant G'(M) parameterizes the dissipative friction process 
that drives the magnetization towards an equilibrium 
state.'^ In the most general case, G(M) is a symmetric 
positive definite matrix that depends on the magnetiza- 
tion direction; however, it is often assumed to be inde- 
pendent of M and proportional to the unit matrix, as- 
sumptions which are valid for isotropic systems. Gilbert 
damping is important in magnetization dynamics. It de- 
termines the magnitudes of the external magnetic field^ 
and the current densitie^ that are required to reorient 
the magnetization direction of a ferromagnet. Therefore, 
a thorough understanding of its properties is essential for 
modeling ferromagnetic systems. 

The main contribution to the Gilbert damping process 
in metalli c ferro magnets is the generation of electron- 
hole pairs.'^J^lHSl model that captures this process was 
developed by Kambersky.'^ In this model, the electrons 
are excited by a time-varying magnetization via electron- 
magnon coupling. If the ferromagnet is in metallic con- 
tact with other materials, the spin-pumping into the ad- 
jacent leads provides an additional contribution to the 
magnetization relaxation.El A general theory that cap- 
tures both of these effects was recently developed.'^ The 
model expresses the G'(M) tensor in terms of the scatter- 
ing matrix S of the ferromagnetic system (m = M/Afg): 
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The expression is evaluated at the Fermi energy. Instead 
of G(M), one often parameterizes the damping by the 
dimcnsionless Gilbert damping parameter a = G/jMs- 
Eq. ([2]) allows studying both the effects of the system 



shape and the disorder dependency of the magnetization 
damping beyond the relaxation time approximation.!^ 

In anisotropic systems, the Gilbert damping is ex- 
pected to be a symmetric tensor with non-vanishing off- 
diagonal terms. We are interested in how this tensor 
structure influences the dynamics of the precessing mag- 
netization in (Ga,Mn)As. Therefore, to briefly discuss 
this issue, let us consider a homogenous ferromagnet in 
which the magnetization direction m = mo -I- Sm pre- 
cesses with a small angle around the equilibrium direc- 
tion mo that points along the external magnetic field 
Hext-'^ For clarity, we neglect the anisotropy in the 
free energy and choose the coordinate system such that 
mo = (0 1) and Sm = {rux my 0). For the lowest order 
of Gilbert damping, the LLG equation can be rewritten 
as: rh — —7m x Hoxt +7m x (a [Hoxt x m]), where a [...] 
is the dimcnsionless Gilbert damping tensor that acts on 
the vector H^xt ^ m- Linearizing the LLG equation re- 
sults in the following set of equations for nix and m^: 
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Here, a,-? are the matrix elements of a when the tensor 
is evaluated along the equilibrium magnetization direc- 
tion mo. For the lowest order of Gilbert damping, the 
eigenvalues of ^ are A± — ±i^Hcxt — iH^xtCt, and the 
eigenvectors describe a precessing magnetization with a 
characteristic life time r = (a7i?oxt)~^- The effective 
damping coefficient a 
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The value of a is generally anisotropic and depends on 
the static magnetization direction mo. The magnetiza- 
tion damping is accessible via ferromagnetic resonance 
(FMR) experiments by measuring the linear relation- 
ship between the FMR line width and the precession fre- 
quency. This linear relationship is proportional to 
and thus depends linearly on a. Therefore, an FMR ex- 
periment can be used to determine the effective damping 
coefficient a. In contrast, the off-diagonal terms, ai^* 

and a^J , do not contribute to the lowest order in the 
damping and are difficult to probe experimentally. 

In this paper, we use Eq. ([2| to study the anisotropy 
and disorder dependency of the Gilbert damping in the 
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ferromagnetic semiconductor (Ga,Mn)As. Damping co- 
efficients of this material in the range of a ~ 0.004 — 0.04 
for annealed samples have been reported.'^^l'SlThe damp- 
ing is anisotropically dependent on the magnetization 
direction. ^-'^ ■^'^ The few previous calculations of the 
Gilbert damping consta nt in this material have indicated 
that a - 0.003 - 0.04.EIilHni These theoretical works 
have included the effects of disorder phenomenologically, 
for instance, by applying the relaxation time approxima- 
tion. In contrast, Eq. ^ allows for studying the disor- 
der effects fully and non-perturbatively for the first time. 
In agreement with Ref. 1151 we show that spin-pumping 
from the localized d-electrons to the itinerant holes dom- 
inates the damping process in the clean limit. In the 
diffusive regime, the breathing Fermi-surface effect is bal- 
anced by effects of the interband transitions, which cause 
the damping to saturate. In determining the anisotropy 
of the Gilbert damping tensor, we find that the shape of 
the sample is typically more important than the effects 
of the strain and the cubic symmetry in the GaAs crys- 
tal.'^ This shape anisotropy of the Gilbert damping in 
(Ga,Mn)As has not been reported before and provides 
a new direction for engineering the magnetization relax- 
ation. 



II. MODEL 

The kinetic-exchange effective Hamiltonian approach 
gives a reasonably good description of the electronic 
properties of (Ga,Mn)As.'i^^ The model assumes that the 
electronic states near the Fermi energy have the character 
of the host material GaAs and that the spins of the itin- 
erant quasiparticles interact with the localized magnetic 
Mn impurities (with spin 5/2) via the isotropic Heisen- 
berg exchange interaction. If the s-d exchange interac- 
tion is modele d by a mean field, the effective Hamiltonian 



takes the formP^lgfll 
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where -ffnoics is the k • p Kohn-Luttinger Hamiltonian de- 
scribing the valence band structure of GaAs and h(r) • s 
is a mean field description of the s-d exchange interac- 
tion between the itinerant holes and the local magnetic 
impurities (s is the spin operator). The exchange field 
h is antiparallel to the magnetization direction m. The 
explicit form of i?Hoios that is needed for realistic model- 
ing of the band structure of GaAs depends on the doping 
level of the system. Higher doping levels often require 
an eight-band model, but a six- or four-band model may 
be sufficient for lower doping levels. In the four-band 
model, the Hamiltonian is projected onto the subspace 
spanned by the four 3/2 spin states at the top of the 
GaAs valence band. The six-band model also includes 
the spin-orbit split-off bands with spin 1/2. The spin- 
orbit splitting of the spin 3/2 and 1/2 states in GaAs 
is 341 meV.i^ We consider a system with a Fermi level 
of 77 meV when measured from the lowest subband. In 



this limit, the following four-band model gives a sufficient 
description: 
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Here, p is the momentum operator, Ji are the spin 3/2 
matrice^^U and 71 , 72 and 73 are the Kohn-Luttinger pa- 
rameters. V(r) = J2i V^i^i^ — R-i) is the impurity scat- 
tering potential, where is the position of the impurity 
i and Vi are the scattering strengths of the impurities^^ 
that are randomly and uniformly distributed in the in- 
terval [— Vb/2, Vo/2]. i^strain is a strain Hamiltonian and 
arises because the (Ga,Mn)As system is grown on top 
of a substrate (such as GaAs) ^ The two first terms in 
Eq. ^ have spherical symmetry, and the term propor- 
tional to 73 — 72 represents the effects of the cubic sym- 
metry of the GaAs crystal. Both this cubic symmetry 
terrrpSJ and the strain Hamiltonian^'* are small compared 
to the spherical portion of the Hamiltonian. A numeri- 
cal calculation shows that they give a correction to the 
Gilbert damping on the order of 10%. However, the un- 
certainty of the numerical results, due to issues such as 
the sample-to-sample disorder fluctuations, is also about 
10%; therefore, we cannot conclude how these terms in- 
fluence the anisotropy of the Gilbert damping. Instead, 
we demonstrate that the shape of the system is the dom- 
inant factor influencing the anisotropy of the damping. 
Therefore, we disregard the strain Hamiltonian Hstrain 
and the term proportional to 73 — 72 in our investigation 
of the Gilbert damping. 




(Ga,Mn)As 




FIG. 1: We consider a (Ga,Mn)As system attached along the 
[010] direction to infinite balHstic GaAs leads. The scattering 
matrix is calculated for the (Ga,Mn)As layer and one lattice 
point into each of the leads. The magnetization is assumed to 
be homogenous. In this paper, we denote the [100] direction 
as the X-axis, the [010] direction as the y-axis and the [001] 
direction as the z-axis 

We consider a discrete (Ga,Mn)As system with 
transverse dimensions € {17,19,21} nm, Lz € 
{11, 15, 17} nm and Ly = 50 nm and connected to infinite 
ballistic GaAs leads, as illustrated in Fig. [T] The leads 
are modeled as being identical to the (Ga,Mn)As system, 
except for the magnetization and disorder. The lattice 
constant is 1 nm, which is much less than the Fermi wave- 
length Af ~ 10 nm. The Fermi energy is 0.077 eV when 
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measured from the lowest subband edge. The Kohn- 
Luttinger parameters are 71 — 7.0 and 72 = 73 = 2.5, 
implying that we apply the spherical approximation for 
the Luttinger Hamiltonian, as mentioned above.l^ We 
use |h| — 0.032 eV for the exchange-field strength. To 
estimate a typical saturation value of the magnetization, 
we use Ms = 10\j\hx/aQ^p^^ with x = 0.05 as the doping 
level and acaAs as the lattice constant for GaAs.l^^ 

The mean free path / for the impurity strength Vq is 
calculated by fitting the average transmission probability 
T = (G) /Gsh to T{Ly) = l/{l + Ly),[23 where Gsh is the 
Sharvin conductance and (G) is the conductance for a 
system of length Ly. 

The scattering matrix is calculated numerically us- 
ing a stable transfer matrix method. The disorder ef- 
fects are fully and non-perturbatively included by the 
ensemble average (a) — '^nLi'^n/Nj, where Nj is the 
number of different impurity configurations. All the 
coefficients are averaged until an uncertainty S (a) = 

^J[{a^) - (af) /Ni of less than 10% is achieved. The 

vertex corrections are exactly included in the scattering 
formalism. 



III. RESULTS AND DISCUSSION 

Without disorder, the Hamiltonian describing our sys- 
tem is rotationally symmetric around the axis parallel 
to h. Let us briefly discuss how this influences the 
particular form of the Gilbert damping tensor For 
clarity, we choose the coordinate axis such that the ex- 
change field points along the z-axis. In this case, the 
Hamiltonian is invariant under all rotations Rz around 
the z-axis. This symmetry requires the energy dissipa- 
tion E oc m^arh (rh"'" is the transposed of rh) of the 
magnetic system to be invariant under the coordinate 
transformations r — R^r (i.e., (rii )'^a rh — m-^drh 
where m ^ Rzin and a is the Gilbert damping ten- 
sor in the rotated coordinate system). Because a only 
depends on the direction of m, which is unchanged 
under the coordinate transformation R^, a — a and 
R^aRz = a. Thus, a and Rz have common eigenvectors 
{|±) = ± i |y))/\/2, and the spectral decompo- 
sition of a is a = q;+ |+) (+| + a_ |— ) (— I + |z) 
Representing the damping tensor in the {\x) , \y) , |z)} 
coordinate basis yields a = dtxx — oiyy ~ («+ + a-)/2, 
otzz — oizi oiyx — OLxy — 0, and a+ — a_. The last 
equality results from real tensor coefficients. However, 
azz cannot be determined uniquely from the energy dis- 
sipation formula E cx rh-^cim because rh is perpendicu- 
lar to the z-axis. Therefore, azz has no physical signifi- 
cance and the energy dissipation is governed by the single 
parameter a. For an infinite system, this damping pa- 
rameter does not depend on the specific direction of the 
magnetization, i.e., it is isotropic because the symmetry 
of the Hamiltonian is not directly linked to the crystallo- 
graphic axes of the underlying crystal lattice (when the 
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FIG. 2: (a) The dimensionless Gilbert damping parameter 
a as a function of the magnetization direction for a system 
where Lx ~ 17 nm, Ly = 50 nm and Lz = 17 nm. (b) The 
dimensionless Gilbert damping parameter a as a function of 
the magnetization direction for a system where = 21 nm, 
Ly = 50 nm and = 11 nm. Here, 9 and (j) are the polar and 
azimuth angles, respectively, that describe the local magne- 
tization direction m — (sin cos smOsm(f>, cos 9). In both 
plots, the mean free path is ? ~ 22 nm 



cubic symmetry term in Eq. (|6| is disregarded). For a 
finite system, the shape of the system induces anisotropy 
in the magnetization damping. This effect is illustrated 
in Fig. [2] which plots the effective damping in Eq. Q as 
a function of the magnetization directions for different 
system shapes. When the cross-section of the conductor 
is deformed from a regular shape to the shape of a thin- 
ner system, the anisotropy of the damping changes. The 
magnetization damping varies from a minimum value of 
around 0.004 to a maximum value of 0.006, e.g., the 
anisotropy is around 50%. The relaxation process is 
largest along the axis where the ballistic leads are con- 
nected, i.e., the y-axis. This shape anisotropy is about 
four- to five-times stronger than the anisotropy induced 
by the strain and the cubic symmetry terms in the Hamil- 
tonian ([6]), which give corrections of about 10 percent. 
For larger systems, we expect this shape effect to be- 
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come less dominant. In these systems, the anisotropy 
of the bulk damping parameter, which is induced by the 
anisotropic terms in the Hamiltonian, should play a more 
significant role. The determination of the system size 
when the strain and cubic anisotropy become comparable 
to the shape anisotropy effects is beyond the scope of this 
paper because the system size is restricted by the com- 
puting time. However, this question should be possible to 
investigate experimentally by measuring the anisotropy 
of the Gilbert damping as a function of the film thickness. 




FIG. 3: The effective dimensionless Gilbert damping Q as a 
function of the disorder. Here, I is the mean free path and 
Ly is the length of the ferromagnetic system in the transport 
direction, cimin and amax are the minimum and maximum val- 
ues of the anisotropic Gilbert damping parameter and Qmcan is 
the effective damping parameter averaged over all the magne- 
tization directions. The system dimensions are = 19 nm, 
Ly — 50 nm and L^ — 15 nm. 

We next investigate how the magnetization relaxation 
process depends on the disorder. Ref. [15] derives an ex- 
pression that relates the Gilbert damping parameter to 
the spin- flip rate T2 of the system: a oc r2(l -I- (T2)^)^^. 
In the low spin- flip rate regime, this expression scales 
with T2 sls a (X T2^, while the damping parameter is 
proportional to a oc T2 in the opposite limit . As ex- 
plained in Ref. 1151 the low spin-flip regime is dominated 
by the spin-pumping process in which angular momen- 
tum is transferred to the itinerant particles; the trans- 
ferred spin is then relaxed with a rate proportional to 
Tg"^. This process appears inside the ferromagnet itself, 
i.e., the spin is transferred from the magnetic system to 
the itinerant particles in the ferromagnet, which are then 
relaxed within the ferromagnet. Therefore, this relax- 
ation mechanism is a bulk process and should not be 
confused with the spin-pumping interface effect across 
the normal metal | ferromagnet interfaces reported in 
Ref. [71 In (Ga, Mn)As, this bulk process corresponds to 
spin-pumping from the d-electrons of the magnetic Mn 
impurities to the itinerant spin 3/2 holes in the valence 
band of the host compound GaAs. The transfer of spin 
to the holes is then relaxed by the impurity scattering 



within the ferromagnet. By contrast, the opposite limit 
is dominated by the breathing Fermi-surface mechanism. 
In this mechanism, the spins of the itinerant particles are 
not able to follow the local magnetization direction adia- 
batically and lag behind with a delay time of 72. In our 
system, which has a large spin-orbit coupling in the band 
structure, we expect the spin-flip rate to be proportional 
to the mean free path {I cx T2).'^° The effective dimen- 
sionless Gilbert damping Q is plotted as a function of 
disorder in Fig. |3j The damping (amoan) partly shows 
the same behavior as that reported in Ref. (121 For clean 
systems (i.e., those with a low spin- flip rate regime), the 
damping increases with disorder. In such a regime, the 
transfer of angular momentum to the spin 3/2 holes is the 
dominant damping process, i.e., the bulk spin-pumping 
process dominates, amcan starts to decrease for smaller 
mean free paths, implying that the main contribution 
to the damping comes from the breathing Fermi-surface 
process. Refs. ,11.16.17. have reported that the Gilbert 
damping may start to increase as a function of disorder in 
dirtier samples. The interband transitions become more 
important with decreasing quasi-particle life times and 
start to dominate the intraband transitions (The intra- 
band transitions give rise to the breathing Fermi-surface 
effect). We do not observe an increasing behavior in the 
more diffusive regime, but we find that the damping sat- 
urates at a value of around 0.0046 (See Fig. [3|. In this 
regime, we believe that the breathing Fermi-surface effect 
is balanced by the interband transitions. The damping 
does not vanish in the limit 1/1 = due to scattering 
at the interface between the GaAs and (Ga,Mn)As lay- 
ers in addition to spin-pumping into the adjacent leads 
(an interface spin-pumping effect, as explained above). 
Fig. [3] shows that the shape anisotropy of the damp- 
ing is reduced by disorder because the difference between 
the maximum (c^max) and minimum (amin) values of the 
damping parameter decrease with disorder. We antici- 
pate this result because disorder increases the bulk damp- 
ing effect, which is expected to be isotropic for an infinite 
system. 



IV. SUMMARY 

In this paper, we studied the magnetization damping 
in the ferromagnetic semiconductor (Ga,Mn)As. The 
Gilbert damping was calculated numerically using a 
recently developed scattering matrix theory of mag- 
netization dissipation.'Sl We conducted a detailed non- 
perturbative study of the effects of disorder and an inves- 
tigation of the damping anisotropy induced by the shape 
of the sample. 

Our analysis showed that the damping process is 
mainly governed by three relaxation mechanisms. In the 
clean limit with little disorder, we found that the magne- 
tization dissipation is dominated by spin-pumping from 
the d-electrons to the itinerant holes. For shorter mean 
free paths, the breathing Fermi-surface effect starts to 
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dominate, which causes the damping to decrease. In 
the diffusive regime, the breathing Fermi-surface effect 
is balanced by the interband transitions and the effective 
damping parameter saturates at a value on the order of 
0.005. 

For the small samples considered in this study, we 
found that the shape of the system was typically more 
important than the anisotropic terms in the Hamiltonian 
for the directional dependency of the damping parame- 
ter. This shape anisotropy has not been reported before 



and offers a new way of manipulating the magnetization 
damping. 



V. ACKNOWLEDGMENTS 

This work was partially supported by the European 
Union FP7 Grant No. 251759 "MACALO". 



^ For a review, see D. C. Ralph and M. Stiles, J. Magn. 
Magn. Mater. 320, 1190 (2008), and reference therein. 

^ B. Heinrich, D. Fraitova, and V. Kambersky, Phys. Status 
Solidi 23, 501 (1967); V. Kambersky, Can. J. Phys. 48, 
2906 (1970); V. Korenman and R.E. Prange, Phys. Rev. 
B 6, 2769 (1972); V.S. Lutovinov and M.Y. Reizer, Zh. 
Eksp. Teor. Fiz. 77, 707 (1979) [Sov. Phys. JETP 50, 355 
(1979)]; V.L. Safonov and H.N. Bertram, Phys. Rev. B 61, 
R14893 (2000); J. Kunes and V. Kambersky, Phys. Rev. B 
65, 212411 (2002); V. Kambersky Phys. Rev. B 76, 134416 
(2007). 

^ T.L. Gilbert, IEEE Trans. Magn. 40, 3443 (2004). 

* J.A.C. Bland and B. Heinrich, Ultrathin Magnetic Struc- 
tures III Fundamentals of Nanomagnetism (Springer Ver- 
lag, Heidelberg, 2004). 

5 V. Kambersky Czech. J. Phys. B 26, 1366 (1976). 

® K. Gilmore, Y.U. Idzerda, and M.D. Stiles, Phys. Rev. 
Lett. 99, 027204 (2007). 

Y. Tserkovnyak, A. Brataas and G.E.W. Bauer, Phys. Rev. 
Lett. 88, 117601 (2002). 
® A. Brataas, Y. Tserkovnyak, and G.E.W. Bauer, Phys. 
Rev. Lett. 101, 037207 (2008); Phys. Rev. B 84, 054416 
(2011). 

^ A. A. Starikov, P. J. Kelly, A. Brataas, Y. Tserkovnyak, 
and G. E. W. Bauer, Phys. Rev. Lett. 105, 236601 
(2010); Y. Liu, Z. Yuan, A. A. Starikov, and P. J. Kelly 
arXiv:1102.5305 

° A similar analysis is presented in J. Seib, D. Steiauf, and 
M. Fahnle, Physical Review B 79, 092418 (2009). 

^ J. Sinova, T. Jungwirth, X. Liu, Y. Sasaki, J.K. Furdyna, 
W.A. Atkinson, and A.H. MacDonald, Phys. Rev. B 69, 
085209 (2004). 

^ Y.H. Matsuda, A. Oiwa, K. Tanaka, and H. Munekata, 

Physica B 376-377, 668 (2006). 
^ A. Wirthmann et ai, Appl. Phys. Lett. 92, 232106 (2008). 
-'^ Kh. Khazen et al, Phys. Rev. B 78, 195210 (2008). 
^ Y. Tserkovnyak, G.A. Fiete, and B.L Halperin, Appl. 

Phys. Lett. 84, 5234 (2004). 
^ I. Garate and A.H. MacDonald, Phys. Rev. B 79, 064403 



(2009). 

L Garate and A.H. MacDonald, Phys. Rev. B 79, 064404 
(2009). 

This shape anisotropy in the Gilbert damping should not 
be confused with the shape anisotropy (in the anisotropy 
field) caused by surface dipoles in non-spherical systems. 
T. Jungwirth, J. Sinova, J. Masek, J. Kucera, and A. H. 
MacDonald, Rev. Mod. Phys. 78, 809864 (2006). 
M. Abolfath, T. Jungwirth, J. Brum, and A. H. MacDon- 
ald, Phys. Rev. B 63, 054418 (2001). 

P.Y. Yu and M. Cardona, Fundamentals of Semicon- 
ductors: Physics and Materials Properties, 3rd Edition 
(Springer Verlag, Berlin, 2005). 

Note that in Eq. ([6| the spin operator s (in the p-d ex- 
change term) is represented in the basis consisting of the 
four spin 3/2 states (s — J/3). The factor 1/3 is absorbed 
in the exchange field h. 

In the discrete version of Eq. ([6|, as used in the numerical 
calculation, we have one impurity at each lattice site. 
A. Chernyshov, M. Overby, X. Liu, J.K. Furdyna, Y. 
Lyanda-Geller, and L.P. Rokhinson, Nature Physics 5, 656 
(2009). 

A. Baldereschi and N.O. Lipari, Phys. Rev. B 8, 2697 
(1973). 

The prefactor of 10 comes from 4 Ga atoms per unit cell 
times spin 5/2 per substitutional Mn, which are assumed to 
be fully polarized. The reduction of the net magnetization 
due to the interstitial Mn ions and p holes are disregarded 
in our estimate. 

S. Datta, Electronic Transport in Mesoscopic Systems 

(Cambridge University Press, Cambridge, England, 1995). 

T. Usuki, M. Saito, M. Takatsu, R. A. Kiehl, and N. 

Yokoyama, Phys. Rev. B 52, 8244 (1995). 

D. Steiauf and M. Fahnle, Physical Review B 72, 064450 

(2005). 

The spin relaxation time of holes in GaAs is on the scale of 
the momentum relaxation time. See D.J. Hilton and C.L. 
Tang, Phys. Rev. Lett. 89, 146601 (2002), and references 
therein. 



